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Abstract 

A list decoding algorithm for matrix-product codes is provided when 
Ci , . . . , C s are nested linear codes and A is a non-singular by columns 
matrix. We estimate the probability of getting more than one codeword as 
output when the constituent codes are Reed- Solomon codes. We extend 
this list decoding algorithm for matrix-product codes with polynomial 
units, which are quasi-cyclic codes. Furthermore, it allows us to consider 
unique decoding for matrix-product codes with polynomial units. 

1 Introduction 

Matrix-product codes, [C\ ■ ■ ■ C s \ - A, arc a generalization of several classic codes 
constructions of codes from old ones [21 [H]. For instance, they extend the 
(u\u + w)-construction. An algorithm for unique decoding when the codes are 
nested, C\ D ■■■ D C s , and A has a certain property, called non-singular by 
columns, was provided in [3]. The algorithm decodes up to half of the minimum 
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distance, assuming that we have a decoding algorithm for d that decodes up 
to half of its minimum distance, for every i. 

List decoding was introduced by Elias [5J and Wozencraft (T5J. The list 
decoder is a relaxation over unique decoding that allows the decoder to produce 
a list of codewords as answers. It can uniquely decode beyond half of the 
minimum distance in some cases or to produces a list of codewords. 

In 1997, Sudan presented a polynomial time algorithm for decoding low 
rate Reed-Solomon codes beyond the classical | bound. Later [7J, Guruswami 
and Sudan provided a significantly improved version of list decoder which can 
correct codes of any rates. Recently, Lee and O'Sullivan provide a list decoding 
algorithm based on the computation of a Grobner basis of a module [13] and 
Beelen and Brander provide an algorithm that has linear complexity in the code 
length [I]. 

In this paper we consider a list decoding algorithm for matrix-product codes 
which is an extension of the algorithm in [5]. The algorithm in [5] assumes a 
known decoding algorithm for every constituent code Ci that decodes up to half 
of its minimum distance, for this algorithm, we assume that the decoding algo- 
rithm is a list-decoding algorithm. Moreover, it is also required that C\,...,C a 
are nested and A is non-singular by columns. The extension is natural, but, we 
believe, it is non-trivial task since we had to modify the algorithm to deal with 
lists of codewords, compute the error bound r and prove the correctness of the 
algorithm, among others. 

Matrix-Product codes are generalized concatenated codes [H 2] , which have 
an efficient decoding algorithm [BJ. However, this algorithm cannot be success- 
fully applied if the matrix A is small, as it is in practice for Matrix-product 
codes (see Remark [374]) . 

The probability of getting more than one codeword as output of a list decod- 
ing algorithm for Reed-Solomon codes was bounded in [T~4]. In section^ we use 
this computation to estimate an upper bound of the probability of getting more 
than one codeword as output, when Ci, . . . , C s are Reed-Solomon codes. The 
algorithm in section[3]can become computationally intense, an optimal situation 
arises considering s = I = 2. 

In section [5] we extend the algorithm in section [3] for matrix-product codes 
with polynomial units which are quasi-cyclic codes. Quasi-cyclic codes be- 
came important after it was shown that some codes in this class meet a modified 
Gilbert- Varshamov bound [TU] , however there are no general fast algorithms for 
decoding them. In [9], many of these codes with good parameters were ob- 
tained. Using list decoding of matrix-product codes with polynomial units we 
can uniquely decode these codes up to the half of the minimum distance. 

2 Matrix-Product Codes 

A matrix-product code is a construction of a code from old ones. 

Definition 2.1. Let C±, . . . , C s C F™ be linear codes of length m and a matrix 
A = (o>i,j) G A4(¥ q ,s x I), with s < I. The matrix-product code C = 
[C\ ■ ■ ■ C s ] ■ A is the set of all matrix-products [ci ■ • • c s ) ■ A where Ci G C; is an 
m X 1 column vector Cj = (ci^, . . . , c m ,i) T for i = 1, . . . , s. Therefore, a typical 
codeword c is 
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/ ci.ioii H hci-o.i ••■ ciiai.H hci, a a s i 



(1) 



\ c mj iai,i + h c m)S a Sj i ••• c m ,iai,z + • • • + c m ,sa s ,i 

Clearly the i-th column of any codeword is an element of the form X^=i a j',i c j € 
F™, therefore reading the entries of the m x Z-matrix above in column- major 
order, the codewords can be viewed as vectors of length ml, 



(2) 



From the above construction it follows that a generator matrix of G is of the 
form: 

^ ai^Gi ai^Gi ■ ■ ■ ai >a Gi ■ ■ • a\jG\ \ 

02,1^2 02.2G2 • • • Cl2,sG2 ' • • CL2JG2 

\ fls,iG s a s ,2G s ■ ■ ■ a StS G s ■ ■ ■ a s jG s ) 
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where Gj is a generator matrix of Gj, i = 1, . . . , s. Moreover, if d is a [m, fcj, dj] 
code then one has that [Gi • • • G s ] • A is a linear code over ¥ q with length and 
dimension k = k\ + • • • + k s if the matrix A has full rank and k < k% + • • • + k s 
otherwise. 

Let us denote by Ri = (cii,i, • ■ • , <Xj,;) the element of consisting of the i-th 
row of A, for i = 1, . . . , s. We denote by Di the minimum distance of the code 
Gi^ generated by (Ri, . . . ,Ri) in F^. In [TS] the following lower bound for the 
minimum distance of the matrix-product code G is obtained, 

d(C) > min{d 1 D 1 ,d 2 D 2 , d s D s }, (3) 

where di is the minimum distance of G^. If Gi, . . . , C s are nested codes, G\ D 
■ ■ ■ D C s , the previous bound is sharp [8] . 

In [2], the following condition for the matrix A is introduced. 

Definition 2.2. [2] Let A be a s x Z matrix and A t be the matrix consisting of 
the first t rows of A. For 1 < ji < ■ ■ ■ < jt < I, we denote by A(ji, . . . ,j t ) the 
t x t matrix consisting of the columns ji , . . . , j t of A t . 

A matrix A is non-singular by columns if A(j\, . . . ,jt) is non-singular for 
each 1 < t < s and 1 < j\ < • • • < j t < I. In particular, a non-singular by 
columns matrix A has full rank. 

Moreover, if A is non-singular by columns and C\ D ■■■ D C s , we have 
d(C) = min{Wi, (I - l)d 2 , ...,(/- s + l)d s } 0. 

In [5] were presented a decoding algorithm for the matrix-product code G = 
[Gi • • • C s ] ■ A C F™ ; , with A non-singular by columns and C\ D ■ ■ ■ D C s , 
assuming that we have a decoding algorithm for Ci, for i = l,...,s. The 
algorithm in [8] decodes up to half of the minimum distance. In next section we 
provide a list decoding algorithm for such codes, assuming that we have a list 
decoding algorithm for Gj, i = 1, . . . , s. 
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3 List Decoding Algorithm for matrix-product 
codes 



Let CcFJ and r > 1. For r G F™, a list decoding algorithm with error bound 
r provides a list with all codewords in C that differ from r in at most r places. 

If r < LhtJ> ii; wiU result int0 

unique decoding. 

We present a list decoding algorithm for a class of matrix-product codes, it 
is an extension of [5J Algorithm 1]. Namely, we consider s nested linear codes 
Ci,...,C s C F™ and a non-singular by columns matrix A G A4(¥ q ,s x I), 
where s < I. We provide a list decoding algorithm for the matrix- product code 
C = [Ci • • ■ C s ] ■ A C F™ , assuming that we have a list decoding algorithm 
LDCi for Ci with error bound Tj. In particular, each LDCi answers an empty 
list if there is no codeword in C,; within distance n of the received word. 

Our list algorithm for C decodes up to 

r = minora + (I - 1), (/ - l)r 2 + (I - 2), . . . , (I - s + l)r s + I - s}. (4) 

We first describe the main steps in our decoding algorithm. The algorithm 
is outlined as a whole in procedural form in Algorithm [T] 

Consider the codeword c = Q^ =1 dj,iCj, ■ • ■ , 2j=i a j-i c j)i where cj G Cj, 
for all j. Suppose that c is sent and that we receive p = c + e, where e = 
(ei, e2, . . . , ei) G F™' is an error vector. We denote by pi = J2j=i a j,i°j + e « e ^™ 
the i-th block of p, for i = 1,. . . J. Let . . . ,i s } C {1, . . . , 1} be an ordered 
subset of indices. We now also suppose that e satisfies the extra property that 

wt{e ij ) < Tj for all j G {1, . . . , s}. (5) 

Since C\ D ■■■ D C s , each block Ylj=i a j,i c j of c is a codeword of C\. 
Therefore, we decode the i\-Xh block pi 1 of p using LDCi and we obtain a list 
L\. Since wt^e^) <t\, we have Ylj=i a j-ii c j ^ ^i- ^ n practice we do not know 
which one of the elements in L\ is J2"j=i a j-i c ji therefore we should consider the 
following computations for every element in L\. Assume now that we consider 
Sj=i a j-ii c j S Li, hence we obtain = pi 1 — Y]j—i a j,h c j an d we can eliminate 
ci in every other block (although we do not know c\) in the following way: we 
consider a new vector p 2 ) G F™ z with components 

Pi =Pt —(Pii - e 4l ) = V af -Cj + e i; for i ^ ti, 

U 

where = Oj,i — ~^- a j,in and p 2 ^ = pi x — e< 1 . Since A is a non-singular 
by columns matrix, the elements of the first row of A are non-zero, and so the 
denominator oi,,^ is non-zero. 

Since C% D ■ ■ ■ D C s , we notice that the i-th block of p 2 ' is a codeword 
of Ci plus the error block a, for i G {1, . . . , s} \ {ii}- We now decode the 
i 2 -th block pf 2 = J2 S j=2 a j,i 2 c i + e *2 0I " P 2 ^ usm g LDC2 and we obtain a list L 2 . 
Since w(ei 2 ) < r 2 , we have X)j= 2 a fio c j ^ ^2- 111 practice we do not know again 

which one of the elements in L 2 is 5Zj=2 a ^\ 2 c j> therefore we should consider the 
following computations for every element in L 2 . Assume now that we consider 
^2 S j = 2 a ^\ 2 c j ^ ^ 2 > hence we obtain ei 2 and, as before, we can eliminate c 2 in 
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every other block (although we do not know C2) as follows: we consider a new 
vector p 3 ) G F™ J with components 

e i2 ) = ^ a % c i + e " for * 7^1, i 2 , 
i=3 

, 3) 2) a 2) 2) 3) 2) , 3) 2) 

where a A = - -vf-a^, = v t [ and p l2 = Pl > - e h . 

2,^2 

Notice that the i-th block of p 3 ' is a codeword of C3 plus the error block e^, 
for i G {1, ...,s}\{ix,i 2 }. 

Then we iterate this process, defining p fe ' for k = 3, . . . , s, and decoding 
the ifc-th block using LDCk- In this way, we obtain the error blocks e^, and 
the corresponding codeword blocks J2j=i a j,i c ji f° r * ■ • • j *s}- The vector 

(J~)j—i ajj x Cj, . . . , Y)j—i Q-j,i s Cj) formed from these s decoded blocks is equal to 
the product [c\ ■ ■ ■ c s \ ■ A(i\, . . . , i s ), where A(i\, . . . , i s ) is the s x s-submatrix of 
A consisting of the columns i\, . . . , i s . Since this matrix is full rank, we can now 
easily compute c\, . . . , c s by inverting A(ii, . . . ,i s ) or solving the corresponding 
linear system. Finally we recover the remaining I — s codeword blocks "for free" 
(i.e. no decoding procedure is involved for these blocks) by simply recomputing 
the entire codeword c = [c\ ■ ■ ■ c s ] ■ A = (J2j=i a i-\ c h •••■> £j=i a j,i c j)i since we 
know the Cj's and the matrix A. 

For each elimination step in the above procedure, it is necessary that 7^ 
0, for each k = 2, . . . , s, to avoid zero division. We claim that this follows from 
the non-singular by columns property of A, exactly in the same way as in [8]. 
Let A 1] = A. The matrix A k) = (a^) G M(¥ q , s X I), k = 2, . . . , s, is obtained 
recursively from A k ~ 1 * > by performing the following I — (k — 1) elementary column 
operations: 

fc-i) 

column* (A fe) ) = column, (A k - 1] ) - "*~ M column^ (A^), 

a k-i,i k -i 

for each i £ {ii, . . . , ik-i}- These operations introduce l — (k — l) additional zero 
elements in the k — 1-th row of A k ^ at each iteration. Hence the submatrix of A k>> 
given by the first k rows and the i\, . . . , ik columns, is a triangular matrix (in this 
case, a column permutation of a lower triangular matrix) whose determinant is 
0*1^ • ■ ■ a k\ k ■ Since A is non-singular by columns, this submatrix is non-singular. 

It follows that the determinant is non-zero, and therefore af} ■ ^ 0. 

The procedure described above will generate a list that includes the sent 
word, if for a given choice of indices . . . , i s } C {1, . . . , each error block 
satisfies wt(ey) < Tj, for all j — l,...,s. The output's procedure may not 
include the sent word if wt(eij) > Tj for some j. 

In the previous description, we have only shown the computations for one 
of the different choices that the list decoder algorithms d, for i = 1, . . . , s, may 
give us. However, if #L\ > 1 then we should consider a different word p 2 ^ for 
every £ G L\. One can see how this tree is created for every element in L in 
line 8 of Algorithm [T] If a decoder LDCj outputs an empty list, for all possible 
choices in Lj—±, then we consider another ordered subset of indices, and start 
the procedure again. 



3) 2) u 2,i , 2) 

Pi =Pi ~ -^TKPii - 

a 2,* 2 
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We now prove that for every error vector e with wt(e) < r there exists a 
good set of indices {i\, . . . ,i s } C {1, . . . , 1} satisfying condition ([5]). We should 
repeat the procedure described above, with every ordered subset of indices and 
collect all the decoded words, in order to be sure that the "good" set of indices 
is considered. 

Theorem 3.1. Let C be the matrix-product code [C\ ■ ■ ■ C s }-A, where G\ D • • • D 
C s and A is a non-singular by columns matrix. Let e = (e\, ei, ■ ■ ■ > ej) € F™ z be 
an error vector with wt(e) < r (see Then there exists an ordered subset 

{ii,...,i s } C{1,...,1} satisfying wt(e i:j ) < Tj, for all j G {l,...,s}. 

Proof. We claim that there exists i\ such that wt(ei 1 ) < t\. Suppose that there 
is no i\ G {1, ...,/} with wt(ei 1 ) < n, that is, wt(ei) > t\ + 1, for alH = 
This implies that 

wt(e) = wt(ei) + • • • + wt{e{) > It\ + I > r 

which contradicts our assumption. 

Let us assume that the property holds for a subset {ii, . . . , ij-x} C {1, . . . , 1} 
of size j — 1 < s. We now prove it holds for a subset of size j. Suppose 
that there is no ij with wt(ei j ) < Tj, that is, wt(ei) > Tj + 1, for all i £ 
{1, . . . , 1} \ . . . , ij-i}- This implies that 

3-1 i 

k=l k=j 
3-1 

> wt ^ ) + ( l -j + i ) T J + ( l -j + i ) 

k=l 

> (l-j + l)Tj + (l-j + l)>T 



which contradicts our assumption and the result holds. □ 

Summarizing, we can now formulate our decoding algorithm for C = [C± ■ ■ ■ C s ] ■ 
A C F™' , where C\ D ■ ■ ■ D C s and A is a non-singular by columns matrix, in 
procedural form in Algorithm [T] 

Corollary 3.2. If wt(e) < r then c is in the list given as output of Algorithm 
[H Hence, the algorithm described in this section is a list decoding algorithm 
with error bound r, i.e. L = {c 6 C \ wt(p — c) < t}. 

Proof. By 13.11 there exists an ordered subset {ii,...,i 8 } C {l,...,l} satis- 
fying wt(ei j ) < Tj, for all j G {1, . . . , s}. Therefore, Cj. G Lj and c = 
(Ylj=i a jA c j> • • • i a j-i c j) e ^- Furthermore, all the words at distance r 

from the received word are included in the output list as well. □ 

Example 3.3. Consider the matrix-product codes with matrix A of the form 




and C\ D C2 Rccd-Solomon Codes over Fig with parameters [15, 10, 6] and 
[15,4,12], respectively. Therefore, the code C — [C1C2] • A has parameters 
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Algorithm 1 List decoding algorithm for C = [C\ ■ ■ ■ C s ] • A 

Input: Received word p = c + e with c € C and wt(e) < r. C\ D ■ ■ ■ D C s 

nested codes and A a non-singular by columns matrix. Decoder LDCi for 

code Cj, i = 1, . . . , s. 
Output: List of all codewords that differ from p in at most r places. 

p' = p; A' = A; Dec = {}; 
for {h,...,i 3 } C {1, ...,1} do 
p = p'; A = A'; U' = {p}; 
for j = 1, . . . , s do 
U = U'- U' = {}; 
for u in U do 
L = LDC 3 { Ul] )- 
for i in L do 

tmp= (0, ...,0) eF m ; 
for k = j + 1, . . . ,s do 
tmp ik = u lk - ^f±£; 

end for 

U' = U' U {trap}- 
end for 
end for 

if U' = {} then 

Break the loop and consider another i\, . . . ,i s in line 2; 
end if 

for k = j + 1, . . . , s do 

column^, (A) = column^ (A) — ^^columni . (A); 

end for 
end for 
for u in U' do 

Obtain (ci, . . . , c s ) from , . . . , Ui s ; 
p = [ci • • • c s ] ■ A; (see (TTJ) and ©) 
if wt(p — p') < r then 

Dec = DecU {p}; 
end if 
end for 
end for 



[30,14,12]. We have error bounds t\ = 3, t 2 = 7, for C\ and C*2 respectively, 
using the list decoding algorithm in [I] or in jTS] with multiplicity w = 4 (see 
next section for further details). Therefore, the error bound for Algorithm [T] is 
t = 7. Note that the error correction capability of C with [5] is only t = 5. 

Let c = (0,0) be the sent word and p = (a 2 x + ax 5 + a 5 x e + a u x 13 ,a 5 x 2 + 
a x e + a s x 10 ) the received word, i.e. wt(e) = 7. 

• We consider the ordered set of indices {1, 2}. So, we decode p\ = a 2 x + 
ax 5 + a 5 x 6 + a 14 x 13 with the list decoding algorithm for C\\ we obtain 
PY = a 2 x + ax 5 + a 5 x 6 + a u x 7 + a w x 13 + a 5 x u . 

Then we compute p 2 } = P2 — Pi = ct 2 x + a 5 x 2 + ax 5 + a 5 x 6 + ax 7 + 
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a 8 x 10 + a 10 a; 13 + a 5 x 14 and decode it with the list decoding algorithm for 
C2. However, we get an empty list as output and we do not consider any 
codeword for the final list. 

• We consider now the ordered set of indices {2, 1}. Therefore, we decode 
P2 = a 5 x 2 + a 7 x 6 + a 8 x 10 with the list decoding algorithm for C\. Wc 
obtain as output p 2 ' = 0. Thus, we compute p\ — p 2 = Pi and we decode 
it with the list decoding algorithm for C2. We have obtain 0. Therefore, 
we deduce that the sent codeword is (0, 0). 

Remark 3.4. Matrix-product codes are generalized concatenated codes 
There is an efficient decoding algorithm for generalized concatenated codes [3], 
the cascaded decoding of multilevel concatenations. In [6], Guruswami and 
Rudra generalize this algorithm to a list decoding algorithm for generalized 
concatenated codes. 

The latter are defined as follows: consider s outer codes, say C 3 out over W q «j 
with parameters (N, Kj, Dj) for j = 0, . . . , s — 1. Let C out = C® ut x ■ ■ • x C s a ~l = 
{(c°, . . . , c s_1 ) I cP 6 C 3 out ,j = 0, . . . , s — 1}, understanding c 3 G C 3 out as a row 
vector. Thus, a typical element c £ C ou t is a s x N matrix, we denote by Ck the 
fc-th column of c, for k = 0, . . . , N — 1. 

One also considers a inner code Ci n over ¥ q and a one-to-one map ij) from 
¥ q a x ■ ■ • x F g « s _i to d n that maps (ioj ■ • ■ , i s -i) 6 F g <« x • • • x F 5 «,_i to a 
codeword ^(io, ■ • ■ , i s -i) hi Cin- 

A generalized concatenated code V of order s is the set 

V = {^(c ), ■ • ■ , V(civ-i) I (co, • • ■ , cjy-i) G B}. 

Guruswami and Rudras's algorithm |6] works as follows (we do not attempt 
to write the algorithm, just to describe its main idea): let R G A4(ao + • • ■ + 
a s _i x N,¥ q ) be the received word. For j = 1, ... , s — 1 consider the code 
Cf n generated by all the rows of the generator matrix in Ci n except the first 

ao H hOj-i, say G 3 in , where C? n = Cj n . The algorithm assumes the existence 

of list decodable algorithms for C{ n and list recovery algorithm for C J out . It is 
also assumed that the list decoding algorithm for C\ n returns a list of messages 
while the list recovery algorithm for C J out returns a list of codewords. 

In the first round one applies, for i = 0, . . . , N — 1, a list decoding algorithm 
of C°„ to 'ipici) obtaining a list Sf. For i = 0, . . . , JV — 1 and for each element 
in Sf, one can recover a message Cj and project it into the first component, 
obtaining a list if. Then it applies list recovery algorithm of C® ut to {Tf}i 
obtaining a list L$ which in particular contains c . 

In the second round we proceed as follows, for each c G Lq, c G .M(ao x 
iV, F g ), consider the matrix c G A^(ao + • • • + a a _i x N, ¥ q ) which is the matrix 
of zeroes with c in the first ao rows. One may subtract R = R — (G°„)*c, i.e., 
if c = c° we are cancelling it from the received word. Thus we have a new 
received word R that should be decode with the generalized concatenated code 
with outer codes C\ ut , . . . , C s Q ~l and inner code C\ n . Since the number of outer 
codes has dropped by one, repeating this process s-times one can successfully 
list decode the original generalized concatenated code. 

Let Ci n be a linear code over ¥ q with parameters [n, k, d] and generator ma- 
trix A, and let C 3 out be a linear code over F g with parameters [N, Kj,Dj] for j = 
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0, . . . , k — 1. Consider ip the encoding linear map of Ci„, i.e., i[){io, • ■ • , ik-i) = 
(io, . . . , Then the generalized concatenated code V is equal to the matrix 

product code [Cl ut , . . . , C^ ut ] ■ A. Since a matrix-product code [Ci, . . . , C s ] ■ A, 
with A non-singular by columns, is a generalized concatenated code one might 
use the algorithm in [6j for matrix-product codes as well. However, the algo- 
rithm in [5] can not be successfully applied for matrix-product codes because 
the inner code has generator matrix A, that is a small matrix (in practice). 

Theorem 13 . 1 1 guarantees the existence of a good set of indices, i.e., satisfying 
< Tj, for all j £ {1, . . . , s}. Hence, in the worst case, we may have 
to consider s!( ) iterations. However, in average we will consider much fewer 
iterations. Given a fix set of ordered indices S — {ix, . . . ,i s }, we will estimate 
for how many error patterns of weight t, one has a good set of indices. In other 
words, what is the probability that a set of indices {ii, . . . , i s } is a good in the 
worst case, i.e. when r errors occur. 

Proposition 3.5. If r errors occur the probability that a fix set of indices 
{ii, . . . , i s } verifies that wt{e.i- ) < Tj for all j G {1, . . . , s} is: 

ETl Y^ min { r ~ a l. T 2} _mm{T-X]^ = 1 a j' T ^-l} fm\fm\ I m \ / »(<-!+!) \ 

Proof. Let us compute the different error vectors e with weight r that allow us 
to have a good set of indices. If wt(ei 1 ) < tx, then we may have < ax < Ti 
errors in block ix, and there are (™) possibilities for having ai errors in block 
ix- Assuming that ax errors occurred in the block ix, we may have wt(ei 2 ) < T2 
if and only if there are 02 errors in the block ii, where < a-i < min\r — a\,Ti\ , 
since wt(e) = r. Hence, there are (™) different possibilities for having a% errors 
in the second block. Repeating this argument for the first s — 1 blocks, we may 
have T—y~]jZx a>j errors in the I — s + 1 remaining blocks (including block i s ) and 

therefore there are possibilities for the remaining blocks. Overall 

we have 

m \ fm(t - s + 1)\ 
a s -xj \t - J2jZ{ a i ) 

error patterns that make {ix, ■ ■ ■ ,i s } & good set of indices. The result holds 
since there are error vectors of weight r. □ 

In Example 13.31 we have s = I = 2 and t\ = 3. Therefore the probability 
that either {1,2} or {2,1} is a good set of indices is 

( 1 o 5 )( 1 7 5 ) + ( 1 1 5 )(a 5 ) + ( 1 2 5 )( 1 5 i ) + ( 1 3 5 )( 1 4 5 ) 1 

(30) 2 

Finally we consider the complexity of alogithm[TJ 

Theorem 3.6. Let LDC\, . . . ,LDC S be the list decoding algorithms consid- 
ered in Algorithm^ with error bounds tx,---,t s , respectively, and that output 



n rnin{r — a\,T-2} rain{r Y2j=i a j ^ T s- 1 } / \ / \ 

ee- e : : 

O1 =0 a 2 =0 a,_i=0 v l/ x ' 
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a list with size bounded by D\ , . . . , D s . We denote by Ri the complexity of the 
algorithm LDCi . Then, algorithm [7] has complexity 

Proof. In the worst case, one should consider every ordered set of s elements 
within the I possible indices, that is s!( ) ordered sets. For a fix ordered set, we 
run LDCi which yields a list of size at most D\, in the worst case. For each 
element in this list we run LDCi producing a list of size D2, in the worst case. 
Hence we will have D\Di words that should be decoded with LDC3. Repeating 
this process, in the worst case, we will decode D± ■ ■ ■ D s -i words with LDC S . □ 



4 List decoding of Matrix-product codes from 
Reed-Solomon codes with small s 

The previous algorithm can become computationally intensive as the number of 
blocks I, the number of blocks that we may need to decode at each iteration s 
and the error bounds t\,...,t s , increase. Therefore, there are two interesting 
situations: considering few blocks and considering codes and error bounds such 
that there is a small probability of getting a list with more that one element as 
output of the list decoding algorithms LDCi. 

Let us consider that the constituent codes C\ D ■ ■ ■ D C s arc Rccd-Solomon. 
This family of codes is especially interesting in this setting for two reasons, con- 
sider an [m, k, d] Reed-Solomon code, there is an efficient list-decoding algorithm 
[7] for decoding up to r" errors, with multiplicity v G N, which is computed as 
follows 



1 1 



(Y) , (r,-l)(fc-l) 



— 1, where 



and r„ is calculated so that 



(v+l\ 



< 



1 



< 



r v + 1 



In particular, one has the algorithm in |13j with complexity 0(D i vm 2 ) and 
the one in [T] with complexity 0(D A vm log 2 m log log m), where D is the list 
size, v the multiplicity and m the length of the code. Notice that, by fixing v, 
we fix t v and bound the list size D < l v /(k — 1). Thus, one may obtain the 
complexity of algorithm [T] as a function of the multiplicities and the length of 
the constituent Reed-Solomon codes by Theorem 13.61 

Furthermore, one has a bound for the the probability p T v (C) that the output 
of Guruswami-Sudan's algorithm for the code C with error bound t v has more 
than 1 codeword |14j , given that at most t v errors have occurred and assuming 
that all the error patterns have the same probability. It turns out that this 
probability may be very small in practice, for example a Reed-Solomon code 
over F 2 e with parameters [64,20,45] and r 1 = 23, this probability is 10 -25 . 
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We consider codes and error bounds in such a way that p Ti (Cj) is small, 
we abbreviate r^* to Tj. With the notation of the previous section, consider a 
received word p = c + e where e is the error vector with wt(e) < t. Consider 
the ordered set of indices . . . , i s } C {1, . . . , I}, if wt{ei j ) < Tj for every j 
then we say that the set of indices is good (otherwise we say that it is bad). For 
a good set of of indices, the sent word c is in the output list by Theorem 13.11 
Furthermore, we claim that with a high probability the output list just contains 
this word: LDC\(pi 1 ) is going to give as output a list containing X)j=i a j-^ c ji m 
Algorithm[TJ In practice this list will only have one element since the probability 
of getting just one codeword is 1 — p Tl (Ci). Then we eliminate c\ in the block 
%2 and decode it using LDC2- Since w{ei 2 ) < T2 we obtain a list that contains 
J2 S j=2 a fi 2 c j an< ^ with high probability this list has only one codeword. We 
proceed in the same way for the rest of the blocks and with probability 

s 

H(i-ptM)) (6) 

i=l 

we obtain an output with just one codeword for this set of indices. 

Consider now a bad set of indices {ii,...,i s }, that is, there exists j such 
that w(ei 1 ) < n, . . . , w(ei j _ 1 ) < Tj-\, but w{ei j ) < Tj, then the block j will 
not be correctly decoded. Again, with probability (O we will obtain at most 
one codeword p' for this set of indices, we do not know anything about this 
codeword excepting that it is not the sent one. However, we claim that with 
a high probability the codewords obtained with this bad set of indices will be 
discarded in line 26 of Algorithm [I] namely, we claim that wt(p — p') > r with 
at least probability 1 — lp Tl (Ci). 

Lemma 4.1. Let p, p' 6 C = [C\ ■ ■ • C s ] ■ A, with p ^ p' and C\ D ■ • • D C s 

Reed-Solomon codes. For t as in we have 

POt(p-p') <r) <lp T1 {Ci). 
Proof. If wt(p — p') < r then there is j such that 

wt(pj - p'j) < t/1 < t/s < T\. 
One has that P{wt(pi — p'A < n) = p Tl (Ci), since p ^ p' andpi,p- G C\. Thus, 

P(wt(p - p') < r) < Y, P (™t(Pi-p'i) <Ti) = iPn(Ci) 

i=l 

since p ^ p' and wt(p - p') = £\ =1 wt{pi - p' t ). □ 

An optimal situation arises considering s = I = 2 and two Reed-Solomon 
codes C\ Z> C*2 such that for error bounds t\ and T2, respectively, Guruswami- 
Sudan's Algorithm outputs a list of at most 1 element with a high probability 
(if at most t\, T2, respectively, errors have occurred). When I = s = 2, this 
construction gives the same family of codes as the (w, u + -^-construction, for 
instance Reed-Muller codes are obtained in that way. 
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5 Bounded Distance Decoding of Quasi-Cyclic 
Codes 



Let C\, . . . , C a C F™ be cyclic codes of length m and A = (aij) an s x Z-matrix, 
with s < I, whose entries are units in the ring F 9 [x]/(x m — 1) or zero. A unit in 
F 9 [x]/(x m — 1) is a polynomial of degree lower than m whose greatest common 
divisor with x m — 1 is 1. The so-called matrix-product code with polynomial 
units is the set C = \C\ ■ ■ ■ C s ] ■ A of all matrix-products [c\ • ■ • c s ] ■ A where 
Ci G Ci C F g [a;]/(a;™ 1 — 1) for i = 1, . , . , s. These codes were introduced in [5]. 

We consider always a special set of matrices A to be defined below with full- 
rank over ¥ q [x]/(x m — 1). Let C, with parameters [ra, ki,di], then the matrix- 
product code with polynomial units C = [C\ ■ ■ ■ C s ] ■ A has length Im and 
dimension k = k\ + • • • + k s . 

Let Ri = (a^i, . . . , a iy i) be the element of (¥ q [x]/(x m — l)) 1 consisting of the 
i-th row of A, where i = 1, . . . , s. Let Cr^ be the ¥ q [x]/(x m — l)-submodule of 
(¥ q [x]/(x m — I)) 1 generated by R\, . . . , Ri. In other words, is a linear code 
over a ring, and we denote by Di the minimum Hamming weight of the words 
of C^, Di = min{wt(x) \ x S C^}. In [3] the following bound on the minimum 
distance was obtained 

d(C) >d*= mm{d 1 D u d 2 D 2 , d s D s }. (7) 

One of the differences between matrix-product codes and matrix-product 
codes with polynomial units is that the lower bound d* is not sharp for the 
latter class of codes. 

The minimum distance can actually be much larger than d* and several codes 
with very good parameters were obtained in this way in [5]. We will provide 
a bounded distance decoding algorithm for these codes, using the list-decoding 
Algorithm [TJ 

Matrix-product codes with polynomial units are quasi-cyclic codes |12j of 
length ml. Although this family provides codes with very good parameters 
there are no general fast algorithms for decoding them. The algorithm in [8 
for matrix-product codes may be used for these family of codes under certain 
hypothesis, but in that case it only corrects up to [ d ^ J . 

We remark that the units of a ring form a multiplicative group, however 
they do not form an additive group. That is, if f,g £ F g [x]/(a; m — 1) are units, 
then fg is a unit but / + g or / — g are not a unit in general. This phenomena 
will impose further restrictions for the list-decoding algorithm since we cannot 
divide by a non-unit. 

Definition 5.1. Let A be a s x I matrix, whose entries are units in the ring 
¥ q [x]/(x m — 1) or zero. Let A t be the matrix consisting of the first t rows of A. 
For 1 < ji < ■ ■ ■ < jt < I, we denote by A(ji, ■ ■ ■ ,jt) the t x t matrix consisting 
of the columns jx, . . . , j t of A t . 

A matrix A is unit by columns if the determinant of A(ji, . . . ,j t ) is a unit 
in ¥g[x]/(x m - 1) for each 1 < t < s and 1 < j± < ■ ■ ■ < j t < I. 

In particular, a unit by column matrix is a non-singular by columns matrix. 
Let C be a matrix-product code with polynomial units i.e. C = [C\ ■ ■ ■ C s ] ■ A, 
where Let C\ D • • • D C s and A is a unit by columns matrix, in particular the 
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elements of the first row of A are non-zero. 

With the notation of section [3J consider a received word p = c + e where e 
is the error vector with wt(e) < t. 

For s = 1, the definitions of non-singular by column- and unit by column - 
matrix are the same. Namely, we can use Algorithm [T] without any modifica- 
tions: for {ii} C {1, . . . , Z} a good set of indices, we decode the block pi 1 with 
LDCi because the cyclic codes generated by / and by fu, with / | x m — 1 and 
gcd(u,x m — 1) = 1, are the same code. Then we divide by ai^ to recover c\ 
(in line 24), we can consider the inverse of ai^ since the entries of A are units. 
Actually, this algorithm for s = 1 is the list decoding version of the algorithm 
in PT| for 1-generator 1-level quasi-cyclic codes. 

For s > 2, for each elimination step in Algorithm [TJ we are dividing by aj^ j 
(in lines 11, 20), we claim that this can be performed because ajj j is a unit. Let 
A k ^ denote the matrix obtained recursively from A by performing the following 
I — (k — 1) elementary column operations (see section 13]): 

fc-i) 

(1 

column 4 (A fe) ) = column, (A fc " 1) ) - fc *~ M column tfc _ 1 (A k ~ 1) ), 

a k-l,i k -i 

for each i ^ {ii, . . . , These operations introduce Z — (fc — 1) additional zero 

elements in the k — 1-th row of A k ^ at each iteration. Hence the submatrix of A k *> 
given by the first k rows and the i\, . . . , ik columns, is a triangular matrix (in this 
case, a column permutation of a lower triangular matrix) whose determinant is 
fih l • ■ ■ ol! ■ . Since A is unit by columns, this minor is a unit. Hence, at ■ is a 
unit, since the units form a multiplicative group. 

Thus, we have a list-decoding algorithm with error bound r as in Fur- 
thermore, we can use it also for unique decoding up to the capacity of the code 
if T=L ^j. 

Theorem 5.2. Consider a matrix-product code with polynomial units C — 
[C\ ■ ■ ■ C s ] ■ A, where C\ D • ■ • D C s and A is a unit by columns matrix. Let 
t = \_ d ^ c ^ 1 j , then the list decoding Algorithm\7\is a unique decoding algorithm 
for C. 

Proof. We have seen above that Algorithm [1] can be successfully applied in this 
setting. Hence, this algorithm is a list decoding algorithm, by Corollary 13.21 
In particular, if wt(e) < r the sent word is in the output list. Moreover, since 
T = L 4^j there is no other codeword at distance r from the received word 
and the result holds. □ 

Example 5.3. Let s = 1, Z = 2, and let C\ be the Reed-Solomon code with 
parameters [15, 8, 8] and generator polynomial / = x 7 + a 6 x 6 + a 13 a- 5 + a 12 x 4 + 
ax 3 + a 10 x 2 + a n x + a 13 , where a is a primitive root Fi6. Let C = [C\] ■ A, 
where A = [1, x A + a 5 x 3 + ax 2 + a x + a 14 ], with a € Fi6 a primitive element. 
One has that C is a quasi-cyclic code with parameters [30, 8, 19]. 

The error correction capability of C is t = 9. However, with the algorithm in 
[5] we can only decode up to 7 errors, since d* = 16. Considering a list decoding 
algorithm with multiplicity 2 for C\, we have an error bound n = 4. Hence the 
error correction capability of Algorithm [T] is t = 2ti + 1 = 9 and we have that 
it is a unique decoding algorithm for C . 
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Example 5.4. Let s = 1, I = 2, and let C\ be the Reed-Solomon code with 
parameters [15, 5, 11] and generator polynomial / = x 10 + a 2 x 9 + a 3 x 8 + a 9 x 7 + 
a 6 x 6 + a 14 .T 5 + a 2 x 4 + ax 3 + a 6 x 2 + ax + a 10 , where a is a primitive element 
in Fi 6 . Let C = [Ci] • A, where A = [1, x 3 + a 3 x 2 + a 14 x + a 9 ]. One has that C 
is a quasi-cyclic code with parameters [30, 5, 24], which is the best known code 
in ®. 

The error correction capability of C is t = 11, however, with the algorithm 
in [5] we can only decode up to 10 errors, since d* — 22. Considering a list 
decoding algorithm for C\ with multiplicity v = 1, we have error bound ti = 5 
and Algorithm [T] decodes up to the half of the minimum distance since its 
correction capability is r = 2t± + 1 = 11. However, considering a list decoding 
algorithm for C\ with multiplicity v = 8, we have error bound n = 7. Hence, 
the error bound for Algorithm [1] is r = 2ti + 1 = 15, which is a list decoding 
algorithm for C. 

Example 5.5. Let s = 1 = 2 and consider a matrix A of the form 



with g a unit in F2[cc]/(x™ — 1). One has that A is a unit by column matrix. 

Consider C\ D Ci Reed-Solomon Codes over Fi6 with parameters [15, 13, 3] 
and [15,8,8], respectively. We consider the unit g — x 5 + a 10 x 3 + a 2 x 2 + a 2 . 
One has that the code C = [C1C2] - A has parameters [30, 21, 7]. Hence, its error 
correction capability is t = 3. Let t± = 1,T2 = 3 be the error bounds for C\ and 
C2, for list decoding algorithms with multiplicity 1. Thus, the error bound for 
Algorithm [1] is t = 3 and we have a unique decoding algorithm for C. 

Note that we may consider unique decoding algorithms for C\ and C2 since 
Ti = [ 2 J ' f° r * = 1)2. Hence, in this case, we can reduce the complexity of 
the algorithm by considering unique decoding algorithms. 

6 Conclusion 

In this article we described a list-decoding algorithm for a class of Matrix- 
Product codes, we computed its error bound and complexity. This algorithm 
can become computationally intense, however we show that for small s, I and 
considering Reed-Solomon codes as constituent codes, the algorithm does not 
become computationally intense. Furthermore, we are able to bound the prob- 
ability of getting more than one codeword as output. The main advantage of 
this approach with respect to Reed-Solomon codes is the possibility of con- 
sidering longer codes without increasing the field size and still using the fast 
decoding algorithms [TJ [T3] for the constituent codes. Moreover, we can con- 
sider a bounded distance decoding, that decodes up to half of the minimum 
distance, for Matrix-Product codes with polynomial units, a family with very 
good parameters. 
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